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We investigate the low-temperature thermodynamics of the unitary Fermi gas by introducing a 
model based on the zero-temperature spectra of both bosonic collective modes and fermonic single- 
particle excitations. We calculate the Helmholtz free energy and from it we obtain the entropy, the 
internal energy and the chemical potential as a function of the temperature. By using these quantities 
and the Landau's expression for the superfluid density we determine analytically the superfluid 
fraction, the critical temperature, the first sound velocity and the second sound velocity. We compare 
our analytical results with other theoretical predictions and experimental data of ultracold atoms 
and dilute neutron matter. 
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I. INTRODUCTION 



II. COLLECTIVE AND SINGLE-PARTICLE 
EXCITATIONS 



In a system of fermions the unitary regime is the situa- 
tion in which r e <C n -1 / 3 <C \a\, where n is total number 
density, r e is the effective radius of the interaction poten- 
tial and a is the s-wave scattering length [jj, |2| . Thus the 
system is dilute but the s-wave scattering length a greatly 
exceeds the average interparticle separation n -1 / 3 . It was 
shown experimentally with dilute and ultracold atomic 
vapors that such systems exist and are (meta)stable Jj|. 
It has been suggested that also the dilute neutron mat- 
ter, which is predicted to fill the crust of neutron stars 
0| , is close to the unitary Fermi gas at a certain den- 
sity range [BJ. At low temperature, the thermodynamic 
properties of the superfluid unitary Fermi gas can be ob- 
tained from the spectrum of elementary excitations, as 
done many years ago by Landau with the superfluid 4 He 
This approach has been adopted by Bulgac, Drut 
and Magierski [Io| and also by Nishida [llj to calculate 
the internal energy and the entropy of the unitary Fermi 
gas. It has been also suggested by Haussmann, Punk 
and Zwerger [l2j], who proposed a way to calculate the 
lifetime of fcrmionic excitations at zero temperature. 

In this paper we adopt the Landau approach by 
introducing a thermodynamical model which uses the col- 
lective bosonic excitations of the generalized hydrody- 
namics [lH and the spectrum of fermionic single-particle 
excitations [3, . We calculate the Helmholtz free en- 
ergy of the two-component balanced unitary Fermi gas 
and from it we determine the entropy, the internal en- 
ergy and the chemical potential. In addition, we use the 
Laudau's criterion to derive the superfluid fraction and 
estimate the critical temperature of the system. Finally, 
by using the obtained superfluid fraction and equations of 
state we calculate the first sound and the second sound of 
the unitary gas as a function of the temperature. Our re- 
sults are compared with previous theoretical predictions 
EE-Ull and experimental data 



For any many-body system the weakly excited states, 
the so-called elementary excitations, can be treated as 
a non-interacting gas of excitations @, Q- In general, 
these elementary excitations are the result of collective 
interactions of the particles of the system, and there- 
fore pertain to the system as whole and not to its sep- 
arate particles @, if- For the unitary Fermi gas the 
mean-field extended BCS theory predicts the existence 
of fermionic single-particle elementary excitations char- 
acterized by an energy gap A [l|, Q ■ The inclusion of 
beyond-mean-field effects, namely quantum fluctuations 
of the orderparameter, gives rise to bosonic collective 
excitations [J 0, which are density waves reducing to 
the Bogoliubov-Goldstone-Andcrson mode in the limit 
of small momenta [l3| . 

The detailed properties of these elementary excitations 
strongly depend on the approximations involved in the 
theoretical approach As previously stressed, in this 

paper we extract the details of the zero-temperature el- 
ementary excitations from a density functional approach 
based on Fixed-Node Diffusion Monte Carlo calculation 
[l3| and from recent Path Integral Monte Carlo simula- 
tions dim. 

It is now well-established [l|, [U, [ll| that the ground- 
state energy Eo of the uniform unitary Fermi gas made 
of N atoms in a volume V is given by 
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where Eq is the ground-state internal energy, £ ~ 0.4 is 
a universal parameter [27j and ep = h 2 (3ir 2 n) 2 / 3 /(2m) 
is the Fermi energy with n — N/V the number density 
and N the number of atoms of the uniform system in a 
volume V. 

The exact dispersion relation of elementary (collective 
and single-particle) excitations is not fully known [l|, . 
In Ref. [1 31 ] we have found the dispersion relation of 
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collective elementary excitations as 



e C oi{q) 



c\q 2 



A 

4m 2 



where 
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is the zero-temperature first sound velocity, with vp = 
(?i/m)(37r 2 n) 1 / 3 the Fermi velocity of a noninteracting 
Fermi gas (see Fig. [IJ. Notice that the term with A 
takes into account the increase of kinetic energy due the 
spatial variation of the density [l3|, [2^-[12- Expanding 
the dispersion relation ([3]) for low momenta we get 



oi(q) = ciq + 
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where the linear term is the familiar phonon dispersion 
relation (the so-called Bogoliubov-Goldstone-Anderson 
mode [H, 0) while the cubic correction depends on both 
the sound velocity c\ and the gradient parameter A. Re- 
cently 



have 



Escobcdo, Mannarclli and Manuel 
eluded an additional dispersive term in the phonon spec- 
trum on the basis of the e expansion of the effective 
field theory [34| . Here the dispersive term depends on 
A. In general, a gradient term with A is essential to 
describe accurately the zero-temperature surface effects 
of a trapped system, in particular with a small num- 
ber of atoms, where the Thomas-Fermi (local density, 
i.e. A = 0) approximation fails [l3j]. For the pourposes 
of the present paper fixing £ = 0.42, that is the Monte 
Carlo prediction for a uniform gas of Astrakharchik et al. 
[35j], we find that the best agreement with Monte Carlo 
data is obtained with A = 0.25. 

The collective modes of Eq. ^ are useful to describe 
correctly only the low-energy density oscillations of the 
system. At higher energies one expects the emergence 
of fermionic single-particle excitations starting from the 
threshold above which there is the breaking of Cooper 
pairs [l], OH [HI- At zero temperature these single- 
particle elementary excitations can be written as 
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where £ is a parameter which takes into account the in- 
teraction between fermions (£ ~ 0.9 according to recent 
Monte Carlo results [l4[) with ep the Fermi energy of the 
ideal Fermi gas. Ao is the zero-temperature gap parame- 
ter with 2Ao the minimal energy to break a Cooper pair 
[HQ. The behavior of e sp (p) is shown in Fig. [TJ where we 
plot also (dotted line) the energy thereshold eth = 2Ao 
above which there is pair breaking and the continuum of 
single- particle excitations [3(| . Expanding e sp (p) around 
the minimum momentum p = \/2mfi = C}^ 2 PFi with 
Pf = \J2mep the Fermi momentum of the ideal Fermi 
gas, we find 

e sp {p) = A + ^-{p-p ) 2 , (6) 
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FIG. 1: (Color online). Elementary excitations of the unitary 
Fermi gas: bosonic collective excitations e co ;(p) (dashed line) 
and fermionic single-particle excitations 2e 3p (p) (solid line). 
The collective mode e co i(p) decays in the single-particle con- 
tinuum when there is the breaking of Cooper pairs, namely 
above eth = 2Aq (dotted line). Zero-temperature parameters 
of elementary excitations: £ = 0.42, A = 0.25, £ = 0.9, and 
7 = 0.45. 



where the effective mass mo is given by 

mAo 



to 
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Notice that the gap energy A of the unitary Fermi gas at 
zero-temperature has been calculated with Monte Carlo 
simulations [HI, HH and reads 7 = Aq/c_f — 0.45. 



III. ELEMENTARY EXCITATIONS AND 
THERMODYNAMICS 

As stressed in the introduction and in the previous sec- 
tion, at very low temperature the thermodynamic prop- 
erties of the supcrfluid unitary Fermi gas can be obtained 
from the collective spectrum given by Eq. ^ and con- 
sidering an ideal Bose gas of elementary excitations [(JIH . 
As T increases also the fermionic single-particle excita- 
tions, given by Eq. ([5]) become important. Thus there 
is also the effect of an ideal Fermi gas of single-particle 
excitations. 

The Helmholtz free energy Fq of the uniform ground 
state coincides with the zero-temperature internal energy 
Eq and is given by 

F = E = ^Ne F . (8) 
5 

The free energy F co ; of the collective excitations is instead 
given by (see also 043) 
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while the free energy F sp due to the single-particle exci- 
tations is 
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Here f3 = l/(fcsT) with T the absolute temperature and 
kg is the Boltzmann constant. The total free energy 
F = F + F co i 



F sp reads 
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where <E>(x) is a function of the scaled temperature x 
T/Tp, with Tp = ep/ltB, given by 
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Notice that the discrete summations have been replaced 
b y integrals, € co1 (t]) = ^^{Xif + 4£/3), and e sp (^) = 
yj [rj 1 — C) 2 + 7 2 - We observe that, by using the expan- 
sions (|4]) and (J6j) for the elementary excitations, adopting 
the Maxwcll-Boltzmann distribution for fcrmionic single- 
particles instead of the Fcrmi-Dirac one, and under the 
further assumption that A = 0, this formula becomes ex- 
actly the simple model 
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proposed Bulgac, Drut and Magierski [T^j. We call this 
equation the BDM model. 

From the Hclmholtz free energy F we can immediately 
obtain the chemical potential /i, that is defined as 
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The chemical potential reads 
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where ^''(x) 



and one recovers /^o = £ef in the 



limit of zero-temperature. 

The entropy S is related to the free energy F by the 
formula 



from which we get 
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FIG. 2: (Color online). Thermodynamical quantities of the 
unitary Fermi gas deduced from our model. Zero-temperature 
parameters of elementary excitations: £ = 0.42, A = 0.25, 
C = 0.9, and 7 = 0.45. 



In addition, the internal energy E, given by 

E = F + TS , 
can be written explicitly as 



E = Ne F 
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To conclude this section we observe that the pressure 
P of the unitary Fermi gas is related to the free energy 
F by the simple expression 
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We can then write the pressure as 
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In Fig. [2] we plot various thermodynamical quanti- 
ties obtained with our model, Eq. (|12j). as a function 
of the scaled temperature T/Tp: the scaled free en- 
ergy F/(Nef), the scaled entropy S/(NkB), the scaled 
chemical potential [iftp and the scaled internal energy 
E/(Ne F ). 



A. Gas of dilute and ultracold atoms 

It is interesting to compare our model, given by Eqs. 
(fTTj) and fF2]) , with other theoretical approaches and also 
with the available experimental data. 

In Fig. [3] we report the data of internal energy E (up- 
per panel) and chemical potential /x (lower panel) ob- 
tained by Bulgac, Drut and Magierski [l(| with their 
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FIG. 3: (Color online). Atomic unitary Fermi gas. Up- 
per panel: scaled internal energy E/(Nep) as a function of 
the scaled temperature T/Tf. Lower panel: scaled chemical 
potential E/(Nef) as a function of the scaled t emp erature 
T/Tf- Filled circles: Monte Carlo simulations [1Q]. Open 
squares with error bars: experimental data of Horikoshi et 
al. [2f|. Solid line: our model, i.e. Eq. JlOj) with Eq. 
p2)l . Dashed line: BDM model [l(J, i.e. Eq. gU) with Eq. 
(|13|1 . Zero-temperature parameters of elementary excitations: 
£ = 0.42, A = 0.25, ( = 0.9, and 7 = 0.45. 



Monte Carlo simulations (filled circles) of the atomic uni- 
tary gas. We insert also the very recent experimental 
data of Horikoshi et al. [25| for the unitary Fermi gas 
of 6 Li atoms but extracted from the gas under harmonic 
confinement (open squares with error bars). In the fig- 
ure we include the results of two models: our model (solid 
line), that is given by Eqs. (JTSJ and ([T2]); the BDM model 
(dashed line), that is given by Eqs. (|TTJ)) and (fT5|) . 

The critical temperature T c of the superfluid-normal 
phase transition has been thoretically estimated to be 
around 0.2Tp. In particular, the theoretical estimations 
are: 0.23 [jul 0.225 0, 0.152 [l7j], 0.15 [3, 0.245 [H, 
and 0.248 [21(. Notice that these values are all much 
smaller than the prediction of the mean-field extended 
BCS theory which is T C /T F = 0.50 0, 0, Gj|. Recent 
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FIG. 4: (Color online). Dilute neutron matter at the density 
n = 0.003 fm~ 3 . Scaled internal energy E/(N€f) as a func- 
tion of the scaled temperature T/Tf- Filled circles: Monte 
Carlo simulations [lij. Solid line: our model, i.e. Eqs. (|15p 
and (TS) with Eq. Dashed line: BDM model QlJ, i.e. 

Eqs. (| 15p and (|19l) with Eq. (I13p . Zero-temperature param- 
eters of elementary excitations: £ = 0.46, A = 0.25, £ = 0.82, 
and 7 = 0.29. 



experiments with 40 K [24| and 6 Li [25[ atoms have mea- 
sured the condensate fraction of the unitary Fermi gas 
and both suggest T c jTp = 0.17. Another very recent 
experiment [26[ has deduced T c /Tp = 0.157 from the 
behaviour of the thermodynamic functions. 

Our model is based on zero-temperature elementary 
excitations and its thermodynamical quantities do not 
show a phase transition. Nevertheless, the results shown 
in Fig. [3] strongly suggests that our model works quite 
well in the superfluid regime, but also slightly above the 
critical temperature (T c ~ 0.15) suggested by two the- 
oretical groups [3, [13] ■ We have also verified that the 
term with A in Eq. (|12p plays a marginal role. The 
main difference between our model and the BDM model 
is instead due to the low-momentum expansions of the 
elementary excitations and to the use of the Maxwell- 
Boltzmann distribution instead of the the Fermi-Dirac 



B. Dilute neutron matter- 
Quantum Monte Carlo data of the dilute neutron mat- 
ter close to the unitarity limit have been produced at fi- 
nite temperature by Wlazlowski and Magierski [l~9| . The 
data have been obtained for the uniform neutron matter 
at the density n = 0.003 mi" 3 , where T F ~ 5-10 10 Kelvin 
(for comparison, Tp ~ 10 -7 Kelvin in ultracold atomic 
vapors). In the neutron matter, the effective radius of 
the nutron-neutron interaction potential is r e ~ 2.8 fm 
and the neutron-neutron scattering length is a ~ —18.5 
fm. This means that in the calculations of Wlazlowski 
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and Magierski pj| r e < n 1 / 3 = d = 6.93 fm < \a\. Thus 
this dilute neutron matter is close but not equal [5J to the 
unitarity Fermi gas (r e <C d <C |o|). @. Consequently, 
the zero-temperature parameters of the elementary exci- 
tations, extracted from the spectral weight function 
are slightly different from those of the unitary Fermi gas 
with a ne glig ible effective range: £ ~ 0.46, £ ~ 0.82, and 
7 ~ 0.29 [ljf- 

In Fig. U]we plot the scaled internal energy E/(Ncf) 
versus the scaled temperature T/Tp of the nuclear mat- 
ter obtained by Wlazlowski and Magierski [l9| with their 
Monte Carlo simulations (filled circles with error bars). 
On the basis of the known zero-temperature parameters 
of the elementary excitations we can compare their finite- 
temperature results with our model (solid line) and the 
BDM model (dashed line) . This value is smaller than the 
one of the atomic unitary Fermi gas because the scaled 
energy gap 7 = A/e^ of the neutron matter at n = 0.003 
fm -3 is smaller than the scaled energy gap of the (atomic) 
unitary Fermi gas. Moreover, the estimated critical tem- 
perature for this dilute neutron matter is T c /Tp ~ 0.09 

In agreement with the findings of Fig. [3l also the re- 
sults of Fig. 2] show that our model (solid line) works 
quite well in the entire supcrfluid regime, but also above 



IV. SUPERFLUID FRACTION 

The total number density n of the unitary Fermi gas 
can be written as 




FIG. 5: (Color online). Superfluid fraction n s /n of the uni- 
tary Fermi g ELS £LS £L function of the scaled temperature T/Tp, 
obtained by using Eq. ([27)) with Eq. p8]l (solid line) and Eq. 
(|29[) (dashed line). Parameters of the elementary excitations: 
£ = 0.42, A = 0.25, £ = 0.9, and 7 = 0.45. 



where the normal density n n , C oi is given by 
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Similar results hold for the normal density n„ iSp due to 
single-particle fcrmionic excitations. It is then easy to 
derive the superfluid fraction 



n = n s + n. 



(22) 
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where n s is the supcrfluid density and n n is the normal 
density At zero temperature n n = and n = n s , 
while at finite temperature the normal density n n is fi- 
nite and increases by increasing the temperature. Corre- 
spondingly, the superfluid density n s decreases and be- 
comes equal to zero at a critical temperature T c . The 
normal density is given by 



(23) 



i.e. the sum of the normal density n„. co ; due to col- 
lective excitations and the normal density n„. sp due to 
the single-particle excitations. According to the Lan- 
dau's approach Q, the gas of collective excitations e co /(p) 
which move with drift velocity v has a distribution 
/b(£co/(p) - p • v), with 



/s(e co /(p)) 
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the Bose-Einstein distribution of collective excitations, 
and total linear momentum 



(25) 



where the universal function of the scaled temper- 
ature x = T/Tf is given by 
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where e coi (r;) = ^yf(\rj 2 +4£/3), and e sp (r?) = 

y/ (v 2 ~~ C) 2 + 7 2 - The function S(x) can be approxi- 
mated as 
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by using the expansions §4§ anc | §Q§ for the elementary ex- 
citations, adopting the Maxwcll-Boltzmann distribution 
for fermionic single-particles instead of the Fermi-Dirac 
one, and assuming A = 0. 

In Fig. [5] we plot the superfluid fraction n s /n of the 
unitary Fermi g function of the scaled temperature 

T/T F , obtained by using Eq. (J27J with Eq. $Z5§ (solid 
line) and Eq. (|29l) (dashed line). The figure shows that 
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the supcrfluid fraction becomes zero at T c jTp = 0.34. 
This value clearly overestimates the critical temperature 
with respect to all other beyond-mean-field determina- 
tions [Hi GJ, [H [Hill HJ. Remarkably the approxi- 
mate formula, Eq. (|29|) . is very close to the full one, Eq. 
(f27|) . up to T/T F ~ 0.15. 



V. SOUND PROPAGATION AT FINITE 
TEMPERATURE 

The analysis of the sound propagation in the supcr- 
fluid unitary Fermi gas at finite temperature can be done 
on the basis of the equations of supcrfluid hydrodynam- 
ics [f| Q , where superfluid and normal densities and ve- 
locities depend on space and time. In our problem the 
constitutive equations to be inserted in the equations of 
supcrfluid hydrodynamics are the Eq. (jlip of the entropy 
S and the Eq. (|2"Tj) of the pressure P. 

According to Landau 6, 8] any superfluid system ad- 
mits a density wave, the first sound, where the velocities 
of superfluid and normal compoments are in-phase, and 
the first sound velocity is given by 



Ul 




(30) 



where S = S/N is the entropy per particle. In addi- 
tion, the superfluid system supports a temperature wave, 
called second sound [H, Q , where the velocities of super- 
fluid and normal components are out-of-phase, and the 
second sound velocity reads 
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Notice that first sound and second sound are given by 
Eqs. (|3U|) and (f5T|) in the hypothesis that these two 
modes are decoupled. As stressed by Taylor el al. 18 1 
this hypothesis is fulfilled as long as R/(R + 1) <C 
(u 2 — u 2 )/{Au 2 u 2 ) where R = (c p ~c v )/c v is the Landau- 
Placzek ratio [39( with c p the equilibrium specific heat 
per unit mass at constant pressure and c v the equilib- 
rium specific heat per unit mass at constant density c v . 
This inequality is met also if R is not small due to the 
fact that the speeds of the first and second sound of the 
unitary Fermi gas are never very close (see below). 

By using our expression (|2ip for the pressure P and 
(§n)s,v = (5/3)P/n [IH the finite-temperature first 
sound velocity becomes 
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From this formula and Eq. (|12j) it is immediate to find 
that for T — > one has u\ — > c\ = vp -\/^/3. By using our 



0.6 



0.5 



0.4 



0.3 



u/v F 




0.4 0.6 

T/T 



FIG. 6: (Color online). Dashed line: scaled first sound veloc- 
ity ui/vf of the unitary Fermi gas as a function of the scaled 
temperature T/T c , obtained using Eq. (|32p with £ = 0.42, 
A = 0.25, 7 = 0.5, and ( = 0.9. Solid line: scaled second 
sound velocity uiIvf of the unitary Fermi gas as a function 
of the scaled temperature T/T c , obtained using Eq. (|33|l with 
£ = 0.42, A = 0.25, ( = 0.9, and 7 = 0.45. 



expression (|17[) for the entropy S the finite-temperature 
second sound velocity can be instead written as 



u 2 = v F 



(*y 



(33) 



From this formula, Eq. (|T2|) and Eq. (|27| with Eq. 

it is not difficult to show that for T — > one has 
U2 — > Ci/v3 = VF^fi/i- In Fig. [6]wc plot the first sound 
velocity u\ and second sound velocity 1*2 as a function of 
the scaled temperature T/T c . These quantities are ob- 
tained by using Eqs. (|3"2"|) and (|33[) . The figure shows 
that ui is weakly dependent on the temperature T while 
U2 strongly depends on T beteween T = and T Cl where 
it vanishes because n s = 0. These results are in quali- 
tative agreement with the recent predictions of Taylor et 
al. based on a T-matrix finite-temperature equation of 
state for the unitary Fermi gas [l8l |. 



VI. CONCLUSIONS 

We have described the elementary excitations of the 
unitary Fermi gas as made of collective bosonic excita- 
tions and fermionic single-particle ones. This approach 
has been used many years ago by Landau with the su- 
perfluid 4 He @ but it is also presently adopted to model 
other many-body systems, like atomic nuclei [40| . Wc 
stress that our approximation of non-interacting elemen- 
tary excitations does not take into account the damp- 
ing of collective modes, which becomes very important 
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by increasing the temperature. We have obtained an 
analytical expression for the Helmholtz free energy and 
the superfiuid fraction, showing that they are sound to 
study the thermodynamics of the unitary Fermi system, 
but only well below the calculated critical temperature of 
the superfiuid phase transition. We believe that this ap- 
proach to the low-temperature thermodynamics can be 
extended to the full BCS-BEC crossover of the Fermi gas 



with two equally-populated spin components. In this case 
the model requires the knowledge of zero-temperature el- 
ementary excitations at finite values of the interaction 
strength 1/ (kpap)- 
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